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Abstract. We present a new approach to cosmological perturbations based on 
the theory of Lie groups and their representations. After re-deriving the standard 
covariant formalism from 50(3) considerations, we provide a new expansion of 
the perturbed Friedmann-Lemaitre-Robertson- Walker (FLRW) metric in terms of 
irreducible representations of the Lorentz group. The resulting decomposition splits 
into (scalar, scalar), (scalar, vector) and (vector, vector) terms. These equations 
directly correspond to the standard Lifshitz classification of cosmological perturbations 
using scalar, vector and tensor modes which arise from the irreducible S*0(3) 
representation of the spatial part of the metric. While the Lorentz group basis matches 
the underlying local symmetries of the FLRW spacetime better than the SO(3), 
the new equations do not provide further simplification compared to the standard 
cosmological perturbation theory. We conjecture that this is due to the fact that 
the so(3, 1) ~ su(2) x su(2) Lorentz algebra has no pair of commuting generators 
commuting with any of the translation group generators. 
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1. Introduction 

Linear cosmological perturbation theory plays a fundamental role in our understanding 
of the evolution of small inhomogeneities of the universe scaling from quantum 
fluctuations to large scale structure formation. In an early paper pQ, Lifshitz and 
Khalatnikov laid the groundwork for all subsequent research on the subject. The 
cosmological perturbations were expressed in spatially flat gauge, concentrating all 
fluctuations into the spatial part of the metric. The perturbations were then classified 
according to their transformation properties in the background spacetime as scalar, 
vector and tensor modes, and it was shown that the different modes do not mix in the 
linear theory. This is the celebrated "decomposition theorem". 

Based on the above results, almost twenty years later, Bardeen introduced a gauge 
invariant formalism [2] which has since become the cornerstone of linear cosmological 
perturbation theory. The formalism has been further developed and generalized in 
subsequent works (see e.g. [3111]). 

In the gauge dependent approach of Lifshitz and Khalatnikov pQ, all gauge freedom 
is used to constrain the perturbations into the spatial domain. Due to the isotropy 
of the background FLRW spacetime, the spatial part of the perturbed metric is an 
SO (3) tensor, and as such, it can be expanded into irreducible representations, an 
I = scalar (the trace of the spatial part of the metric) and an / = 2 tensor (a 
traceless symmetric tensor). The latter has five independent components, corresponding 
to the m = —2,-1,0,1,2 eigenvalues of the rotation generator around an arbitrarily 
chosen direction. The classification of [1] corresponds to \m\ of these representations, 
where \m\ = 0, 1,2 belongs to the scalar-, vector- and tensor modes respectively (note 
the slightly different and confusing terminology from group theory). In other words, 
perturbations are classified according to how a mode responds to rotations around a 
chosen axis. 

Interestingly in [2], Bardeen kept the same classification scheme even for gauge 
invariant perturbations. Once the whole metric is considered, it is locally transforming as 
an SO(3, 1), i.e. a Lorentz tensor. The Lorentz algebra describes transformations of the 
time-time and time-space components of the metric under infinitesimal so(3) rotations 
of space, a subalgebra of the full Lorentz algebra. It is easy to show that the time-time 
component of the metric transforms as a scalar with respect to spatial rotations, and 
the time-space part as a vector. As shown by Bardeen, the decomposition theorem 
still holds using the classification of the perturbations according to the m eigenvalue 
(magnetic quantum number in physicist's terminology). 

The above considerations provide a motivation to investigate the mathematical 
theory of cosmological perturbations from a group theoretical point of view. In addition, 
in the literature it is customary to start from the result that the spatial dependence of 
linear perturbations in the Fourier-space can be completely represented by the solutions 
of the generalized Helmholtz equation. The mathematical background of this result is 
that the symmetry group of the solutions of the Helmholtz equation, the Euclidean group 
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(see e.g. [S]), is precisely the group of geometrical symmetries of the spatial part of the 
(homogeneous and isotropic) FLRW metric. Thus the expansion of the spatial part of the 
perturbations into solutions of the Helmholtz equation is essentially a harmonic analysis 
on the Euclidean group as a Lie group. This observation provides further motivation to 
examine closely the group theoretical structure behind perturbation theory. 

To the best of our knowledge, the theory of cosmological perturbations, in general, 
has not been considered to date from the point of view of Lie groups, and the principal 
goal of the present paper is to develop an approach to cosmological perturbations based 
on the theory of Lie groups and their representations. More specifically, we investigate 
the theory using irreducible representations derived from the more relevant Lorentz 
algebra, and not just its so(3) subalgebra as usual. This is motivated by the observation 
that the spatially flat FLRW spacetime (which we only consider in this paper), i.e. the 
zeroth order cosmological solution, is conformally flat, and by the locally Lorentzian 
nature of the metric tensor. Our hope is to elucidate the role that symmetries play in 
the decomposition of the 10 Einstein equations into independent subsets based on the 
decomposition theorem. 

The plan of the paper is as follows. In section [2, we present a quick overview of 
the standard covariant formalism; for further details please refer to [2j [3j H] or, for a 
more informal introduction, [6] . In section [3j first we reproduce the standard covariant 
formalism from 50(3) considerations, and then we expand the perturbations using 
irreducible representations of the full Lorentz algebra. In section HI we present the 
perturbed field equations obtained from the Lorentz decomposition, and in section [5j 
we discuss gauge invariant quantities. Finally, in section [6j we summarize our results 
and draw conclusions. 

2. Standard covariant formalism in a nutshell 

The background spacetime is described by the FLRW metric which we introduce in the 
conformal form 

ds 2 = g^dx^dx" = a 2 (n)(—dr) 2 + r )ijdx l dx s ) , (1) 

where n is the conformal time variable and 7^ is the metric tensor for a 3-space of 
uniform spatial curvature K. Here, and throughout the paper, Greek indices run from 
to 3 while Latin indices run from 1 to 3. 

As mentioned earlier, perturbations in various quantities are classified according 
to how they transform under spatial coordinate transformations in the background 
spacetime as scalar-, vector- and tensor modes. The homogeneity and isotropy of 
the background metric allows the separation of the time dependence and the spatial 
dependence, with the spatial dependence related to solutions of a generalized Helmholtz- 
equation. Scalar-, vector- and tensor harmonics are solutions of the scalar- 

V 2 Q (0) + A; 2 Q (0) =0, 
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vector- 

V 2 QM+k 2 QW = 0, (2) 

and tensor 

V 2 Qg } + eqf = o 

Helmholtz-equations respectively, where V denotes the covariant derivative with respect 
to the spatial metric 7^. 

In a spatially flat (K = 0) universe, the harmonics are essentially plane waves: 

g (0) = exp(ik ■ x) , 

Qi^ = — ^(ei ±ie 2 )iexp(ik • x) , (3) 
13 

Q\f 2) = - y g(ei ± ze 2 )»(ei ± ze 2 ) 3 - exp(zk • x) , 

where x = (xi,X2,Xs) and ex,e 2 are unit vectors spanning the plane transverse to the 
wave vector k = (k\, /c 2 , /C3). The vector modes represent divergence-free (vorticity) 
vectors while the tensor modes are transverse and traceless, and represent gravitational 
waves: 

V l Ql ±1} = , V'Qif ) = , 7 ^lf } = • (4) 

The curl free vectors and the longitudinal components of tensors can be obtained from 
covariant derivatives of the scalar- and vector modes by 



Qif = (AT 2 V;V,. + - 7y )Q (0) , (5) 



2k 

where k — |k|. 

A completely general perturbation of the gravitational field can be written as a 
linear combination of perturbations associated with individual spatial harmonics defined 
above, and no coupling between the different modes. Thus, a general perturbation to 
the FLRW metric can be represented for example as [2] 

300 = -a 2 (l + 2A), 

9oi = ~ a 2 Bi , (6) 



9ij =a\ li] +2H Lllj + 2H [ 



with 



Tij) , 



A = A(x k , r]) — scalar potential; 
Bi = Bi(x k , n) — vector shift; 

Hl = HL(x k ,rj) — scalar perturbation to the spatial curvature; 
Hxij = HTij(x k , rj) — trace free distortion to the spatial metric. 
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For the kth harmonics, the scalar-, vector- and tensor components of the perturbed 
metric become 

A(77,x) =A( V ,k)Q^(pc), 
H L (r,,x) =H L ( V ,k)Q(°\x.), 

Bib,*) = £ BM(r, : k)Q ( r\x), 



(7) 



m=— 1 
2 



fWu,x)= £ 4 m) (^A;)Q^(x). 



\(m) , 



m=-2 



Let us now rewrite the perturbed metric ([6]) in the form 

9ap = 9a/3 + 8g a p . (8) 

It can be shown by simple calculations that according to the presented expansion, the 
first order part of the metric 

5g 00 = - 2a 2 A , 

Sg Qi = - a 2 Bi , (9) 
Sgij = 2a 2 (H L j ij + H Tij ) , 

with the choice of k = (0, 0, k), can be decomposed as 

*<to= |Q (0) X>A| • do) 



where the Oj amplitudes are 

oi = 2a 2 A(r], k) , 



o 2 = 2a 2 H L ( V ,k) 



o 3 = -^=B^( V ,k), o, = a 2 B(°\ V ,k) 



o 5 = -^B^\ V ,k), (11) 



06 



9 



\c?h£%, k) , o 7 = ^H^\ V , k) , o 8 = ^ifffo *) 



2a_ 2 
IT 



and Oj are the following 10, 4x4, symmetric matrices: 
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This decomposition provides a unique expansion of the metric perturbations up to 
constant normalization factors. 

The parametrization of the matter stress-energy tensor of a perfect fluid in the 
standard formalism is given by 

T° = -p- Sp, 

T c i = {p + p){v i -B i ) 1 

T o = ~{p+P) v i, 

T) = (p + Sp)5 l J +pW J , 

where Sp is a scalar density perturbation, Vi is a vector velocity perturbation to the 
spacelike part of the zeroth order four velocity u Q = [a(r]), 0, 0, 0] of the fluid, and Ily is 
a tensor anisotropic stress perturbation. 

In the present paper we prefer to work with lower indices, thus for the perturbation 
of the stress-energy tensor we use 

8T a p = g a ^ST^ - g av b~g UT T Tf) . (13) 
By taking into account the explicit form of T a p as 

T a/ 3 = (p + p)u a U/3 + pg a p, (14) 
with the given parametrization of the metric and the matter perturbations we get 

5T 00 = a 2 (5p + 2pA), 

5T 0i = -a 2 [(p + p)(v t - Bj+pBi], (15) 
STij = a^dpd^ + plUij + 2{H Llij + Hj<ij)]\. 

Now, in a completely analog fashion to the metric perturbations (given in ( FlOl) ). we 
can decompose the perturbations of the stress-energy tensor as 



10 



(16) 



j'=i 



a/3 



where, for the kth harmonics, the scalar-, vector- and tensor components of the 
perturbations become 

5p(r/,x) =dp( V ,k)Q(°\x), 
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and the 2} coefficients are 



a 2 (5p + 2pA) 



7i 
T 7 



v/2 



(-i) 



v/2 

?! 

2\/2 
a 2 p 



(p + p)v^ 1] - pB 



(p + p)v (+1) - pfi (+1) 



T 2 = a 2 (5p + 2pH L ), 

" 2 T (p + p)v<®-pB® 



s/2 



%a 2 p (n(- 2 ) + 2ff 



i 2 ] 



(n(°) + 2H, 



T j > 



T 10 = f-a 2 p (n( +2 ) + 2« 



(+2) 
T 



7 

(17) 



("2)\ 
T j > 



;is) 



Tg = u> ( n (+u + 2 4 +1) 

2v/2 v T 

With this parametrization it is straightforward, although tedious, to derive the 
Einstein's field equations at linear order. The explicit forms of the scalar-, vector- and 
tensor equations are provided for example in the equations (17), (20) and (22) in [6] 
respectively. We will use them as a reference hereafter, and also present them here in 
the spatially flat (K = 0) case for our latter convenience. 

Scalar equations 



a 1 • o 7i — R 

—A — Hl — -H T = AnGa 2 {p + p) 

A - 



' h , . v — B 
,„ + 3-(, + J ,)_ 



a f a\ 2 ad k 2 
2- -2 - +-^" T 
a Va/ ad?/ 3 



d a 
drj a 



k Pi 1 
(H L + — )=4nGa 2 (5p+-6p). 



Vector equations 



kB {±1) - H { t 1] = 16nGa 2 {p + p) 



;(±1) _ £(±1) 
k 



drj a 



Tensor equations 



(fc S (±D _ h^) = -87rGa 2 prt ±l \ 



d 2 ^ a d ^2 
drj 2 a drj 



H^ 2) = 87rGa 2 prt ±2 \ 
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3. Perturbation theory using irreducible representations 

In the present section, we re-derive the perturbed field equations of the standard 
formalism with more explicit reference to group theory. First we consider the expansion 
of the metric perturbations according to irreducible representation of the £0(3) group, 
and show that the corresponding equations are identical with the ones presented in the 
previous section. After this preparation, we consider the expansion according to an 
irreducible representation of S0(3, 1), the full Lorentz group. We explicitly compute 
a (trace-orthogonal) basis for the decomposition of a second grade S0(3, 1) tensor 
according to the su(2) x su(2) decomposition of the Lorentz algebra. We present the 
corresponding linear field equations in the next section. 

3.1. The 50(3) case 

To specify our notation, the Lorentz algebra is defined by the following commutation 
relations 

[Lpu, L pu ] = ig up Lp a — igppL ua — ig va Tpp + igp a T vp . (19) 

Note that we are using "physicist" convention, i.e. Hermitian generators, instead of the 
ant i- Hermit ian (or "mathematician") convention. The spatial part of the algebra can 
be cast in a familiar form 

[Jj, J k ] = ie jM Ji , (20) 

where Jj = \ejk\Lu is the Hodge dual of the spatial generators. This subalgebra is 
equivalent to that of su(2) or so(3). With the usual boost generators Ki = L 0i we have 
the following relations 

[Kj,K k ] = - ie jk iKi , [Jj,K k ] = ie jk iKi . (21) 

For the FLRW solution in a spatially flat universe, the metric tensor is conformally 
flat and its spatial part is conformally Euclidean. Moreover, the spatial part of the 
metric is a locally SO (3) tensor. Therefore it makes sense to expand the spatial part 
of the metric in terms of irreducible representations under 5*0(3). The irreducible 
representations of a 3 x 3, symmetric matrix are the trace (/ = 0, scalar), and a rank-2, 
traceless, symmetric tensor (1 = 2 mode). Here I corresponds to the eigenvalues of the 
Casimir operator (total angular momentum in physics terminology), and for / = 2 there 
are five modes which correspond to the five degrees of freedom in a rank-2, traceless, 
symmetric matrix. These are usually parametrized according to the eigenvalue of the 
rotation generator in the direction of an arbitrarily chosen z-axis. The corresponding 
eigenvalues (magnetic quantum numbers in physics) are m = —2, —1, 0, 1, 2. These V^ m 
tensors are traditionally called "polarization operators", and we quote their explicit 
forms from [7] as given below. 
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o , 

"2/ 

(22) 

\ 

1 . 
l) 

These matrices yield a natural basis on which to expand 3x3 symmetric matrices, 
such as the spatial part of the metric tensor. Careful comparison reveals that the 
previously introduced 2 matrix corresponds to V^o, i-e. the trace is an 50(3) scalar, 
while Oq . . . Oio directly correspond to V<>,2 • • • Vi -2, the 5*0(3) tensor modes with the 
five possible values of m; the three dimensional representations were extended into four 
dimensions with rows and columns of O's. 

The time-time part of the metric tensor g Q0 does not respond to SO (3) rotations, 
i.e. it is an 50(3) scalar. In four dimensions it is represented with the matrix 0\. It 
can be shown (c.f. eq. (l2Tj) ) that under rotations goi is an 50(3) vector. The spherical 
basis for a vector is g 03 , g i ±Wo2- This is represented with the matrices 4 , 5 and 3 
respectively (up to a — % factor introduced for convenience sake). These latter matrices 
contain zeros for their spatial part. 

As it was noted in the introduction, and it also follows from the above argument, 
the standard parametrization of the metric (based on the solution of the Helmholtz 
equations) directly corresponds to an irreducible representation according to 50(3), 
i.e. spatial rotations. The resulting 10 parameters repackage the 10 degrees of freedom 
in a symmetric 4x4 matrix according to irreducible representations of 50(3). 

These matrices are suitable parameters for an expansion motivated by the 50(3) 
invariance of the spatial part of the FLRW solution at fixed time, and the local 50(3) 
invariance of space. The homogeneity (translation invariance) of the universe can 
be further utilized with Fourier transformation, i.e. an expansion according to the 
irreducible representation of the translation generators with eigenvalue k. In particular, 
once a corresponding fc-mode is fixed for the transform, we can choose the orientation 
of our z-axis (global rotation invariance). This will zero out derivatives in the other two 
directions. Conventionally this idea is expressed using " transverse" and " longitudinal " 
modes in real space. 

We thus recovered the perturbation theory, described in section [2} by expanding 
the perturbed metric according to the above representation of the 50(3) group, and 
Fourier transforming after choosing the z-axis such that it is parallel with the k-mode 
in question. This procedure exploited the maximal symmetry group of the Helmholtz 
equation, the Euclidean group, that is also the group of geometrical symmetries of the 
spatial part of the background solution. The corresponding perturbed field equations 
are identical to the ones of the standard covariant formalism. Using computer algebra, 



2.2 
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-i 1 
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V 1 i °7 
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/ -1 \ 
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we explicitly recovered the equations of [B] based on the matrix formulation presented 
above. 

3.2. The 50(3,1) case 

Let us now take into consideration the larger local symmetry of the metric tensor, 
the SO(3, 1) group, for constructing a perturbation theory according to irreducible 
representations. This is also motivated by the conformal flatness of the spatially flat 
FLRW solution. The calculation is exactly analogous to the exposition of the SO (3) 
perturbation theory above. A popular parametrization of the so(3, 1) algebra splits into 
two commuting su(2) algebras. Then for each su(2) we can use the familiar theory of 
irreducible representations. 

To find the irreducible representations, we used a spinor formalism: in each sw(2)'s 
the irreducible representations are found by the standard way of symmetrizing spinor 
expressions. The direct product of the two commuting irreducible representations form 
a representation of so(3, 1). For constructing spinor representations we used Infeld-van 
der Waerden symbols (see e.g. [8]) 

If 1 \ 1 If 1 \ 

V2 [ 1 ) ' ° AB ' ~V2\1 ) ' 

(23) 

1 / 



a AB' 



a 



AB' 








° AB ' ~ V2 



to map g^ u into Qaa'bb 1 > where primed and unprimed indices live in the two commuting 
su(2)'s. The symmetric 4x4 metric tensor is equivalent to a direct product of two 
vector representations for each of the su(2)'s. One can find matrices analogous to Oj 
by symmetrizing in the primed and unprimed indices. After explicit symmetrization we 
have one scalar mode: the 4D trace (S in ( 124)) ). The traceless 4x4 tensor is a direct 
product of two vector modes, each of which can be classified according to the respective 
to eigenvalue, corresponding to the generator according to an arbitrarily chosen z-axis 
in each representation. The pair (to, to'), both take three possible values, —1, 0, 1, giving 
9 possibilities. These 1+9 modes cover the 10 degrees of freedom in 4 x 4 symmetric 
tensors. We present the corresponding S and S m)m i matrices explicitly: 



Si. 



/ 1 








1 \ 
































V i 








1 ) 




( 








\ 
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i 
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-1 







V o 








o / 



Si: 



0,-1 



/ 1 i \ 

10 1 

% % 

V 1 i / 

/ 1 -i \ 

10 1 

-i -i 

V \ -i J 
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$0,0 



5-1,- 



S-i, 



( 1 \ 
10 
10 

V o o o -l y 

( \ 
1 — i 
-i -1 

V o o o o/ 

/ 1 -1 \ 



V -l o o i / 



<So,i 



>5-i,o 



/ 1 i \ 

10 0-1 

i -i 

V -1 -i / 

/ 1 - /' 
10 0-1 
-i i 

V o -l i o y 

/ -1 \ 
10 
10 

V o o o i y 



(24) 



Similar to the SO (3) case, we can decompose the perturbed metric according to 
the above su(2) x su(2) representation as 



8g a p = « 2 



sS + J2, 



i.'Sn 



"m^m >~"m,m' 
m,m'=—l 



(25) 



where 



Sm,m' Sm,m'(-E )j S s(x ), 



and the stress-energy tensor perturbation as 

i 

<5T nn = -a 2 ' 



ST ni = a 2 



5Tij = a 2 



(t + P s)s + j2(t m,m' P s m,in')S m ^ 

m,m'=—l 
1 

m,m' PSm,m' )S r 

m,m'=— 1 
1 

m,m' PSm,m') Sm,m' 



J 00 



(26) 



0i 



m,m'=— 1 



with 



f. 



m,m' 



tm,m'\-E )? ^ t(x ). 



Once we expand Einstein's equation into the parameters corresponding to the above 
matrices, fix a A:-mode for Fourier transforming, and choose the z-axis parallel with this 
mode, we obtain the perturbation equations. 



4. Perturbation equations 

The decomposition of the perturbations according to the S mjm > and S matrices results 
three independent groups of equations. A simple relation maps our decomposition 
to that of the traditional one, namely m + m! — > m, where (m, m!) corresponds to 
the su(2) x su(2), and m corresponds to the traditional so(3) scalar- vector-tensor 
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decomposition. For example (0,1) or (—1,0) correspond to vector modes, etc. In 
addition, our parameter s also corresponds to a scalar mode in scalar-vector-tensor 
terms. 

Accordingly, the first group consists of two independent wave equations for si,i and 
s_i _i as 



s\ ,i + + k 2 s 1A = 167rGa\ tl , 

a 

2d 

s_i _i H s_i _i + k 2 s^i _x = 167rGa 2 t_i _i 

a 



(27) 
(28) 



where, and throughout the paper, an overdot denotes the derivative with respect to 
the conformal time r]. The above equations are then clearly identical with the tensor 
equations of the standard formalism provided in Sec. 2, and thus the si,i and s_i,_i 
amplitudes correspond to gravitational waves while the ti,i and £_i,_i amplitudes 
correspond to anisotropic stress perturbations. 

The second group of equations consists of two decoupled systems, one for si, and 
so,i, and another for s ,-i an d The two systems are identical. 

K s o,\ + si, ) - i{s ,i - si, ) = 1Qn ^ a [t 0>1 + t lj0 - (p + p)(s ,i + si,o)], (29) 
k( s -ifi + s o -i) - i(s_i, - s -i) = [t-i,o + *o,-i - (P + p)(s-i,o + so -i)], (30) 



d 2d 

dr? a 

d 2d 

dr? a 



(ifc(s ,i + si )0 ) + (s ,i - si )0 ) = 167rGa 2 (i ,i - *i,o), 
(ifc(s_i )0 + s -i) + (s-i,o - s _i) = 167rGa 2 (f_i i0 - *o,-i)- 



(31) 
(32) 



Thus, the si,o, so,i, £i,o, ^o,i an d s_i,o, so,-i, £-i,o, ^o,-i degrees of freedom in 
the su(2) x su(2) representation correspond to the vector amplitudes of the so(3) 
representation. 

The remaining group of equations is the following coupled system of the si,-i, so,o, 
s_i,i, s and ti-i, to,o, £_i,i, t amplitudes, 

3d 2 a 

"^2~( s o,o + s i ,-i + s -i,i ~ s ) + - [«o,o + si -i + s-1,1 + 3s - 2«fc(si _i - s_i,i)] (33) 
+ A: 2 (s , + s) = -87rGa 2 (t ,o + + + * + 2ps), 



d 2d 
dn a 



{h-i + - 2s ,o + 2ik(s-i > i - si _i)) - /c 2 (si _i + s-i,! - 2s) 



(34) 



= 167rGa 2 (t 1 ,_ 1 +t_ 1 , 1 -2too), 



-(s - Sl _l - S-1,1 - so,o) - so,o - s 



SuGa 2 
k 



(*i _i - - (p + p)(si,-i - s_i,i)), (35) 



d 2 a\ d 
a 2 a / dry 



+ fc 2 



(so,o + si -l + s-i,! - s) - 4 



a.. . . . 

-(so,o + si -1 + s_i i + s 
a 



+ (36) 



2ifc 



d d \ 
a dn J 



(si,-i - s_i,i) 



16irGa 2 (t- ps). 



This system is clearly equivalent to the scalar equations of the so(3) representation. 
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In general, one can find the following one-to-one correspondence between the 
standard so(3) and the su(2) x su(2) parametrization of the perturbation amplitudes 




— — \ 7> n T i £-1,-1 




iB^tH^ if. (±1) P ILW\ 

iB^ ± 4 ±1} 1 
s o,±i — ' *±i,o 




{p + p )ivW±?—- , (37) 



(o) 



2^2 ' ±i,U 2v/2 V 2 



h l - a + ifl<°) 4 0) i / p n(°) , 

si,-i = 2 3~' = _ 2 V ° P + °P + ~ § (P + Pj™ 

H L -A-iBW 4 0) . l/ r , pn(°) , (0) \ 

= - 2 | L , t-i.1 = ~2 Up + h + ^y- + (p + pW } ) 

Hl~A2 (o) pn(°) gp+jp 
so,o ^ h 3 ' ,0 ~~ ~ 3 2 

A + 3# L ^ 35p + 5p 

s = ^ , £ = 7. • 

2 2 

Despite the property that the Lorentz group matches the underlying local 

symmetries of the FLRW spacetime better than SO(3), no further simplifications arise 

in the field equations compared to the standard decomposition theorem. We could 

not further decouple the above system. The new equations based on the irreducible 

representation of the SO(3, 1) group through an su(2) x su(2) decomposition of the 

Lorentz algebra appear to be equivalent to the SO (3) approach. 

5. Gauge invariant quantities 

One of the key features of Bardeen's formalism is the systematic determination of 
the gauge invariant variables. It is therefore important to show how these quantities 
can be obtained under the su(2) x su(2) representation. In the present section we 
derive the transformation properties of the metric perturbation amplitudes s and s m>m i, 
and compute the corresponding gauge invariant metric quantities. For the matter 
perturbations the computations are completely analogous. 

The most general gauge transformation is induced by the coordinate transformation 

x a = x a + £ a (x f) ), (38) 

where £ a is of the same order as s and s m ^ m >. 

The changes in the metric tensor can be computed from the transformation law 
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and the scale factors in g a p and g a p at linear order are related by 
a(fj) = a(ri) 



1 + V 
a 



For the linear metric perturbations one can obtain the relation 



8g a p{x<) = 5g af3 (x K ) - a 2 (rj) 



2d 
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(40) 



(41) 



where r]^ is the Minkowski metric tensor, and a comma denotes partial derivative. After 
inserting the explicit su(2) x su(2) decomposed forms of 5g a p and 5g a p using (1231) we 
get 



1 

m,m'=—l 



a 8 



sS -\- ^ Sm,m' S m ,m' 
m,m'=— 1 



'^eVap-V^-V^a- (42) 
(X 

aB 



The final form of the change in the amplitudes is then given by 



Sl,l 


= 5l,l, 


(43) 




= 5-1,-1, 


(44) 


Sl,0 


= 5i,o + \ (-i 1 + - z^ 1 - ke) , 


(45) 


5-1,0 


= 5-1,0 + \ (-i 1 - 1? + *ke - ke) , 


(46) 


50,1 


= 5 ,i + \ (-i 1 + ^e 2 + ike + ke) , 


(47) 


S~0,-1 


= 50,-1 + \ (-e 1 - ii 2 - ike + ke) , 


(48) 


Si -1 


= s 1 ,„ 1 +^(e-e+ike-ike), 


(49) 


5-1,1 


= s- 1 , 1 + ^(e+t 3 -ike-ike), 


(50) 


50,0 


= 5 ,o+^(e°+z^ 3 ), 


(51) 


S 




(52) 



The Si,i and s_i_i amplitudes are gauge invariant, just as we expected since 
they correspond to the tensor modes in the standard formalism. Similarly to the 
field equations, the remaining amplitudes decouple into two groups. The amplitudes 
belonging to the vector modes in the standard formalism depend only on and £ 2 , 
while the amplitudes belonging to the scalar modes depend only on £° and £ 3 . 

Consequently, from the first group, one can construct two gauge invariant variables 



= *(5l,0 + 5 ,l) + 



50,1 - 5l, 



#(-1) 



Hs_i, + sq-x) + 



5-1,0 ~ 5q,_i 

k 



(53) 
(54) 



which coincide with Bardeen's gauge invariant vector amplitudes. 
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From the second group one can again construct two gauge invariant quantities 



$1 

$ 2 



S + S0,0 + - 
a 



2i 

lfe (fl - 1 ' 1 



1 , 

Sl-l) + + s i-i 



2^0 ,0; 



s - so,o - «1 -1 - 

1 



2i 
J 



S-l,l - Sl,-l 



,S_1 i - Si _i 



fc 2 



(55) 



(56) 



2s'o,o ~~ s'i i — S-i,i H — (2so,o — s i,-i ~~ 

a 

which, as one can check, are related to Bardeen's scalar variables as 

$i = 2§h and $ 2 = 2$ A . (57) 

As it is well known from the standard so(3) theory (see e.g. [3]), any gauge invariant 
variable that can be constructed from si >0 , S-i,o, so,i and Sq,-i and their time derivatives 
can be written as a linear combination of and and their time derivatives with 
coefficients of arbitrary functions of time. The same is true for the amplitudes s^-i, 
s -i,i> s o,o and s with the variables $x and <& 2 . In addition, corresponding pair of gauge 
invariant variables can be obtained from the perturbations of the stress-energy tensor to 
form the complete, closed set of gauge invariant variables (up to linear transformations). 



6. Summary and Conclusions 



We have developed an approach to linear cosmological perturbation theory based on 
Lie groups and their representations. After a short overview on the standard formalism, 
based on the decomposition by the solutions of the generalized Helmholtz equation, 
we explicitly reproduced the theory from the irreducible representations of the 50(3) 
group, and utilizing the translation symmetry by Fourier transforming. This procedure 
completely exploited the maximal symmetry group of the Helmholtz equation (i.e. the 
Euclidean group), the group of geometrical symmetries of the spatial part of the 
background FLRW solution. 

As an application of our new approach, we expanded the perturbations according 
to irreducible representations of the Lorentz algebra through the popular su(2) x su(2) 
decomposition. This expansion was motivated by the fact that the spatially flat FLRW 
cosmological metric is conformally flat, and that the local symmetry group of the metric 
is the Lorentz group, of which 50(3) is only a subgroup. Our hope was to elucidate the 
role that local symmetries play in the decomposition of the perturbed field equations, 
and check whether it is possible to obtain a simpler system of equations than that of 
the standard formalism. 

As a result we presented a new decomposition by (scalar, scalar), (scalar, vector), 
and (vector, vector) terms. We showed that these can be mapped directly into the 
standard scalar-vector-tensor classification. Furthermore, we showed that the resulting 
equations are simply related to the equations of [6]. The decomposition according to 
the larger Lorentz group did not produce a finer split of the perturbed field equations 
than the standard one based on the subgroup 50(3). Our calculation elucidated the 
explicit role of the local 50(3) and 50(3, 1) invariance as well as translation invariance. 
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In particular, in the so(3) theory, the J3 generator commutes with the generator of 
translations in the z-direction, a property that is crucial for simplifying the equations. 
On the contrary, in the su(2) x su(2) representation, no so(3, 1) generator commutes 
with any of the translation generators. We speculate that this is the reason why the 
larger group, which most naturally matches the local symmetries of the metric tensor 
and the conformal form of the zeroth order cosmological solution, did not provide a 
simpler system of equations than the standard formulation. 

Our calculation provides a blue-print for decomposing tensor quantities using 
group theoretical considerations in different perturbation theories of Einstein's 
field equations. Possible generalizations of the group theoretical formulation for 
cosmological perturbation theory includes different representations of the Lorentz (or 
the conformal Poincare) group, using for example the full conformal group generators 
for decomposition, fully spinorial perturbation theory, or generalization to higher 
dimensions using the corresponding local symmetries of the metric tensor. The presented 
group theoretical approach could also be useful beyond the linear order, or considering 
inhomogeneous cosmological models. Furthermore the formalism may be applied to 
perturbation theory of other, not necessarily cosmological, solutions of the Einstein 
equations with different background symmetries. 
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